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A recently suggested modified BCS (MBCS) model has been studied at finite temperature. We 
show that this approach does not allow the existence of the normal (non-superfluid) phase at any 
finite temperature. Other MBCS predictions such as a negative pairing gap, pairing induced by 
heating in closed-shell nuclei, and "superfluid - super-superfluid" phase transition are discussed 
also. The MBCS model is tested by comparing with exact solutions for the picket fence model. 
Here, severe violation of the internal symmetry of the problem is detected. The MBCS equations 
are found to be inconsistent. The limit of the MBCS applicability has been determined to be far 
below the "superfluid - normal" phase transition of the conventional FT-BCS, where the model 
performs worse than the FT-BCS. 

PACS numbers: 21.60.-n, 24.10.Pa 



I. MOTIVATION 

Interest in nuclear pairing correlations has been in- 
tensified for many reasons in recent years (see, e.g., re- 
views 21 0] ) . Among the many aspects of the problem, a 
thermal behavior of the pairing correlations in nuclei is 
considered also. It is well known that the conventional 
thermal BCS approach produces not very precise results 
when applied to finite many-particle systems like atomic 
nuclei, principally due to particle number fluctuations. 
In order to overcome, at least partially, the shortcomings 
of the BCS approach a new model, named the modi- 
fied BCS (MBCS), was suggested and explored in papers 
H H H U According to the MBCS calculations a 
sharp "superfluid - normal" phase transition, which is 
a distinct feature of the conventional thermal BCS the- 
ory, appears at much higher temperatures and may be 
smeared out completely. In this paper we analyze the 
performance of the MBCS. 



II. INTRODUCTION TO THE MBCS MODEL 

The conventional BCS theory is based on the Bo- 
goliubov transformation from particle creation aj m and 



Hamiltonian: 



annihilation 
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operators to quasiparticle operators 
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3'~ l jm 



i Vj OLj m , 
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where the index jm corresponds to the level of a mean 
field with quantum numbers j = [n,l,j], projection m 
(we consider the spherical case) and energy Ej. Tilde in 
Eq. and below means time reversal operation: ctj m — 

(-iy-™ a ^ m . 

The BCS equations at zero temperature T = are ob- 
tained, e.g., by minimization of the energy of the pairing 
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in the ground state (treated as a quasiparticle vacuum) 
under the condition that the number of particles N in 
the system is conserved on average. Let us consider the 
simplest case of a constant pairing matrix element Gjj' = 
G. Then the BCS equations have the form: 



iV = 2j2 n i v 



A = Gj2 n . 



jUjVj, 



(3) 



where 




ilj = (j + 1/2) and Ej = y/ (ej — A) 2 + A 2 is a quasipar- 
ticle energy. In these equations A is a pairing gap and A 
is a chemical potential or the energy of a Fermi level. 

The nuclear Hamiltonian H' = i? pa ir — A7V can be 
rewritten in terms of the quasiparticles as 

H' = U + Yj b M + J2 °1 C^J + A o)+ H c + Hres , (5) 
3 3 

where 

N 3 = a ]m a 3™, A] = -j== ^ a] m a]^ 

m V 3 m >0 

and (see e.g. Ref. 0) 
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In Eqs. © the terms which renormalize single-particle 
energies (~ Gv?) are omitted. 

When the pairing strength G is weak, the BCS equa- 
tions yield the trivial solution (normal phase): {uj, Vj} — 
{0(1), 1(0)} for all j, i.e. A = 0. Above some critical 
value G C r, a superfluid solution appears as energetically 
preferable. The value of the pairing gap A, which receives 
a positive contribution from all levels, may be considered 
as a measure of how strong pairing is in the system. In- 
deed, the UjVj combination in the second expression of 
Eqs. © indicates how far away the system is from the 
trivial solution. 

In the conventional BCS theory at finite temperature 
(FT-BCS), minimization of the pairing Hamiltonian is 
replaced by the statistical average of the free energy over 
the grand canonical ensemble. The FT-BCS equations 
read 
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A = (7^J<),: 1 2n j )u j v j , 
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where rij are the thermal Fermi-Dirac occupation num- 
bers for the Bogoliubov quasiparticles 



i, = l/(l + cxp(^/T)) 



(8) 



with Uj and Vj coefficients and energies Ej having the 
same form as Eqs. J2J at T = 0. However, now they are 
temperature dependent through the A and A values. 

The MBCS model also starts from the Hamiltonian 
(Eq. (JSJ) and the canonical Bogoliubov transformation 
(Eq. QJ). New ingredients appear at the extension of 
the approach to finite temperatures. 

In brief, a temperature-dependent unitary transforma- 
tion to the Bogoliubov quasiparticles {c*j m , a jm } is ap- 
plied, thus transforming them into new bar-quasiparticles 

^jm } : 



n j a jm 



(9) 



A new ground state |0) is introduced as a vacuum for the 
bar-quasiparticles 



*jm a 3 m 



|0) = 0. 



(10) 



The coefficients in the transformation, Eq. JJjJ, are se- 
lected so that 



(0> 



15) 



(11) 



and it is assumed that the occupation numbers rij for 
the Bogoliubov quasiparticles should have the same form, 
Eq. JSJ, as in the statistical approach. 

Combining Eq. JIJ and Eq. © the particle opera- 
tors {Oj m , cij m } are expressed in terms of the bar- 
quasiparticle {a! m , otjm\ ones 



t 



a 



jm 



u j a jm + v j a jmi 



(12) 



where 



= UjyJl-rij+Vj^/n, 
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Since the expectation value of -ff pa i r at T 7^ in the 
|0) ground state looks similar in terms of Uj and Vj co- 
efficients to the one at zero temperature in terms of Uj 
and Vj , the MBCS equations are written down in analogy 
with Eqs. © as 



N = 2J2^3 

3 

A — G CljUjVj 



(14) 



or in terms of Uj and Vj coefficients and thermal quasi- 
particle occupation numbers rij 
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III. 



THERMAL BEHAVIOR OF THE MBCS 
PAIRING GAP 



In Refs. 0, H 13], applying the MBCS to study the 
thermal behavior of different nuclear quantities, the au- 
thors point out the following distinctive features of the 
new model: a) the pairing gap decreases monotonically as 
temperature increases and does not vanish even at very 
high T; b) the "superfluid - normal" phase transition is 
completely washed out. 

Taking the MBCS equations as they have been sug- 
gested, we analyze the validity of the above results. We 
have repeated the MBCS calculations for neutrons in Ni 
isotopes and for neutrons and protons in 120 Sn in Ref. 
and Ref. 0, respectively. The MBCS Eqs. have been 
solved with an accuracy of 10 -11 . Our code reproduces 
excellently all the results in Ref. Q. As typical examples, 
we use in this presentation the nuclei 76 Ni (quasibound 
calculation, single particle levels in the continuum having 
no width), 84 Ni (resonant-continuum calculation, finite 
width is taken into account for the levels with positive 
energy) and 120 Sn. 

Starting with the thermal behavior of the pairing 
gap. The neutron MBCS pairing gap in 76 Ni is plot- 
ted in Fig. ^a). One notices that it reaches zero at 
T w 2.1 MeV and continues to decrease with the nega- 
tive sign. These calculations were performed as reported 
in Ref. 3 on a truncated single-particle basis assuming 
the N = — 28 inert core. Later it was recommended 
in Ref. p| that the MBCS calculations should be per- 
formed on an entire or as large as possible single-particle 
spectrum. However, though in Ni isotopes this recipe of 
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FIG. 1: The MBCS pairing gap A (upper part (a)) and Vj 
coefficients (iower part (b)) for particie ieveis near the Fermi 
surface in 76 Ni (neutrons). The conventional critical temper- 
ature T c is shown by the vertical arrow and the single particle 
spectrum used is from Ref. 0- 



"entire spectrum" helps to avoid negative values of the 
pairing gap, it does not work in the case of 120 Sn (see 
below). 

In calculations with a wider single-particle spectrum in 
Ni isotopes we have found that the gap starts to contin- 
uously increase above a certain temperature remaining 
always positive. The pairing strength has been renor- 
malized to keep the A(T = 0) value. 

An example of such behavior of the pairing gap is 
presented in Fig. Ufa) for 84 Ni. As can be seen, at 
T « 3.3 MeV some strange discontinuities are apparent, 
a phenomenon which may be defined as a "superfluid 
super-superfluid" (S-SS) phase transition. At this tem- 
perature, the MBCS equations find a new energetically- 
preferable solution. A T-dependence of the total excita- 
tion energy given by 



E* = £(T) — £(0), 



(16) 



where £{T) is calculated from Eq. (45) in Ref. |4j, is 
shown in Fig. Hfc). The chemical potential jumps away 
from the X(T = 0) value at T w 3.3 MeV (see Fig.^b)). 

In short, while at the A = point the behavior of 
physical observables is smooth (see Fig. ^b)), the dis- 
continuities at the S-SS point in excitation energy, pair- 
ing gap and chemical potential take place. Thus, there is 
no phase transition from the superfluid to normal phase 
within the MBCS but instead a phase transition of a new 
type is predicted at finite temperatures. 
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FIG. 2: The MBCS pairing gap A (upper part (a)), chemical 
potential A (middle part (b)), and excitation energy E* (lower 
part (c)) in 84 Ni (neutrons). Single particle spectrum from 
Ref. 3 is extended by including levels from N = — 28 shells. 



The absence of the normal phase in all MBCS calcu- 
lations has motivated us to apply this model to a magic- 
number system of nucleons in which the existence of this 
phase at T = is expected. The nucleus 120 Sn has been 
taken as an example. The MBCS neutron pairing gap 
(solid line in Fig.|3fa)) in this nucleus shows the same be- 
havior as already discussed for 76 Ni even though a rather 
complete single particle basis has been employed in 120 Sn. 
The neutron pairing gap vanishes at T w 5.5 MeV and 
becomes negative at higher temperatures. 

The proton gap in 120 Sn exhibits a rather strange be- 
havior as a function of T. Starting from zero value at 
T = the gap smoothly develops to a value of —0.73 MeV 
at T = 5 MeV (dashed line in Fig-Efa)). A completely 
different result for the MBCS proton pairing gap in this 
nucleus has been reported in Ref. [j| and it is shown as 
the dotted line Acs m Fig-EJ a )- To obtain it, it has been 
suggested that closed-shell systems should be treated dif- 
ferently from open-shell ones, namely, all summations in 
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FIG. 3: The MBCS pairing gap A (upper part (a)) and 
chemical potential A (lower part (b)) for neutrons (solid lines) 
and protons (dashed lines) in 120 Sn. Acs and Acs (dotted 
lines) correspond to closed-shell calculations in Ref. [Ej. See 
text for details. 



the MBCS equations l|15|) should be carried over hole 
levels only. 

In our opinion, the last recipe contradicts the pre- 
vious recipe of "entire spectrum" in the same publica- 
tion. It is also evident that this artificial constraint 
cannot be justified for a heated system. Indeed, both 
subsystems, the former hole (fully occupied) and parti- 
cle (empty) levels, become partially occupied due to the 
heating and there are no physical reasons to ignore the 
particle part of a single-particle spectrum. Moreover, to 
keep the number of nucleons N constant under the above 
constraint, forces the MBCS equations to abnormally 
renormalize the corresponding chemical potential. In 
the above CS-calculations the A quantity runs away from 
K,cs(T = 0) = -10 MeV to K,cs(T = 5) = +175 MeV 
(following the dotted line in Fig. E^b)). In addition, the 
capacity of the system is zero in the CS-approximation, 
i.e. the system is artificially frozen. 

The above MBCS results are natural consequences of 
Eqs. (|15fl . For example, these equations do not have the 



trivial solution {t 



Eqs. (|T3)) it would correspond to 

lj = ^/l — Uj] 



,} = {0(1), 1(0)}. Indeed, from 



particles 
holes 



j,-! — — y /i,,,, Cj = ^/l — n,j 
and contradict the positive dehnition of Uj. 

One may also notice from Eqs. (|13l) that Vj coefficients 
become negative for particle levels above a certain tem- 
perature, with rij increasing, since Va << uj. In which 
case the MBCS pairing gap A receives a positive contri- 



bution from hole levels and a negative contribution from 
particle levels. Thus, the contribution of single-particle 
and single-hole states to a pairing phenomenon appears 
to be essentially different. 

Numeric calculations (see Fig. ^b)) show that two 
terms in the second expression of Eqs. (| 1 51) compensate 
each other around the critical temperature of the con- 
ventional BCS T c pa 0.57 • Aji =0 for particle levels and 
Vj become negative at higher temperatures. The gap A 
may vanish at some temperature but only when a nega- 
tive contribution from particles and positive contribution 
from holes cancel each other (notice the difference with 
the conventional BCS A = 0). However if this happens, 
at higher T the balance appears to be broken and A be- 
comes finite again. It also means that A = has nothing 
in common with the normal phase and that one cannot 
conclude from the absolute value of the pairing gap how 
strong the pairing is in the system. 

The temperature behavior of the pairing gap depends 
on a delicate balance between particle- and hole-parts 
of the single-particle basis used which makes the MBCS 
predictions very doubtful. 



IV. MBCS AND EXACT SOLUTIONS OF THE 
PAIRING HAMILTONIAN 

In this section we compare MBCS predictions with ex- 
act solutions of the pairing Hamiltonian employing the 
Picket Fence Model (PFM) |9( which is widely used as a 
test model for the pairing problem (see, e.g., Ref. [Io|). 
For numeric calculations we have selected N = 10 levels, 
each of two- fold degenerate (for spin up and spin down), 
with the energy difference of 1 MeV and 10 particles dis- 
tributed over the levels. This configuration thus repre- 
sents 5 levels for holes with energies e_i = —0.5 MeV, 
— 1.5 MeV, etc and 5 levels for particles with energies 
Ei = +0.5 MeV, +1.5 MeV, etc. 

The MBCS predictions for the pairing gap, excitation 
energy, and specific heat given by 



d£_ 
df 



are shown in Fig. 01 by the solid lines. Here, also shown 
are the FT-BCS results as dashed lines (explicit expres- 
sions for the quantity £ in both approaches are given 
below). Both the MBCS and FT-BCS can be compared 
to the exact solutions of the pairing Hamiltonian (dot- 
ted lines). It should be noted that the MBCS-PFM 
pairing gap behavior is qualitatively very similar to the 
one in 84 Ni reported in Fig. Efa) and that the "super- 
fluid - super-superfluid" phase transition takes place at 
T « 1.78 MeV. 

One concludes from Fig. 01 that the MBCS does not 
achieve its main goal of improving upon the description 
of heated nuclei in the conventional FT-BCS. Indeed, ex- 
cept for a narrow region around T c , the deviation from 
the exact results is worse in the MBCS case compared to 
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FIG. 4: Pairing gap A (upper part (a)), excitation energy 
E* (middle part (b)), and specific heat (lower part (c)) C v 
as predicted by the MBCS (solid lines) and FT-BCS (dashed 
lines) for the PFM (N = 10 and G = 0.4 MeV). The exact 
results in (b) and (c) are plotted by dotted lines. 



the conventional FT-BCS, which persists even at a very 
low temperature. It is also obvious that the "superfluid - 
super-superfluid" phase transition as well as the overall 
behavior at higher temperatures are artificial effects of 
the MBCS. 

In addition, a more detailed analysis shows that the sit- 
uation with the MBCS is much worse than the disagree- 
ment apparent in Fig.Q] We further present in Fig.|S]the 
spectroscopic factors for two particle and two hole lev- 
els closest to the Fermi surface, within the MBCS (solid 
lines) and FT-BCS (dashed lines) for a comparison with 
the exact results (dotted lines). It is important to keep in 
mind that the pairing Hamiltonian in the PFM possesses 
particle-hole symmetry. For this reason, dotted curves in 
Fig. [5] for hole —i and particle i levels are ideally sym- 
metric about the y = 1 line. The same is true for the 
FT-BCS results. Nothing of this symmetry remains af- 
ter the secondary Bogoliubov transformation of Eq. J^J is 
applied in the MBCS (see solid lines in the same figure). 
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FIG. 5: Spectroscopic factors of two lowest particle and hole 
levels in the PFM (TV = 10 and G = 0.4 MeV) calculated 
within the MBCS (solid lines) and FT-BCS (dashed lines). 
Dotted lines represent the exact results. 



In addition, the description of this physical observable in 
the MBCS is very poor compared to the FT-BCS results. 

Breaking of the particle-hole symmetry in the MBCS 
calculations is even more clearly seen from an analy- 
sis of the MBCS quasiparticle spectrum presented in 
Fig. Elb). In the PFM it should be two- fold degener- 
ate (i.e. Ei = E-i) as, e.g., in the FT-BCS calculation 
in Fig. Eta) because of this symmetry. The chemical po- 
tential A in the MBCS calculations does not stay at zero 
energy (as it should be) but runs away to positive values 
as temperature increases. This explains why the level 
i = 2 appears at a lower energy than the i = 1 level 
above T 1.6 MeV. It is also the origin of the Ei/E-i 
splitting in these calculations. 

We have repeated the MBCS-PFM calculations with 
different values of the pairing strength G > G CT . As 
G increases, the MBCS "superfluid - super-superfluid" 
phase transition takes place at a lower temperature and 
the splitting of i and —i levels becomes stronger. 

Another consequence of the particle-hole symmetry is 
that under any conditions it should be true that Ui = V—i 
and vi = u_i or alternatively 

u\ + u 2 _, = u\ + vl = 1 . (17) 

In Fig. \7\ we demonstrate what happens to this analyt- 
ical identity in the MBCS. The calculations have been 
performed for a different number of levels N of the PFM 
and the strength parameter G has been adjusted to keep 
T c = 0.5 MeV in each calculation. The convergence of 
the results in this T-range (also for A, A, and E*) is 
reached at N « 10. The conclusion from Fig. [7| is that 
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FIG. 6: Quasiparticle spectrum in the PFM (JV = 10 and 
G = 0.4 MeV) within the FT-BCS (upper part (a)) and the 
MBCS (lower part (b)). Particle (hole) levels are plotted by 
thick (thin) lines. 
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FIG. 7: The (u\ + uLt) quantity in the MBCS-PFM calcu- 
lations with different N. T c — 0.5 MeV is the same in all 
calculations. Analytical identity (y=l) is plotted as a dotted 
line. 



the analytical identity of Eq. (|17[) is completely broken 
in MBCS calculations even from very low temperatures. 

Thus, one can see that the MBCS severely violates the 
symmetry between particles and holes which is an essen- 
tial feature of the pairing problem solved for the PFM. 
Accordingly, the MBCS can be applied to this problem 
only at temperatures much lower than T c . 



V. INCONSISTENCIES OF MBCS AND MHFB 
APPROACHES 

In this section the evaluation of the MBCS equations 
is analyzed to verify consistency. To start with we agree 
that the pairing Hamiltonian H' expressed via the MBCS 
variables Uj and Vj has the same form as the BCS Hamil- 
tonian H' of Eq. JSJ) at T — 0. It is also true that the ex- 
pectation value of the pairing Hamiltonian (0|-ff pa ir|0) at 
T 7^ looks similar to (0|iJ pa i r |0) in the case of the BCS 
at T = 0. But then the formal similarity is broken and 
the corresponding energies Ej are attributed to the con- 
ventional Bogoliubov quasiparticles and not to the "bar"- 

quasiparticles. These energies Ej — y (ej — A) 2 + A 2 en- 
ter through new thermal occupation numbers rij and new 
u, ^-coefficients 




(see remarks on page 8 in Ref. 5] just after Eq. (86)). 

In another words, the MBCS procedure yields new 
eigen energies of Bogoliubov quasiparticles while new 
eigen states are now modified quasiparticles. The point 
here is that if one would take the mean of the modi- 
fied quasiparticle energies under Ej , they should coincide 
with the BCS Ej at T = 0. Indeed, the secondary Bo- 
goliubov transformation of Eq. iJSJ is a unitary one and as 
such, cannot change the eigenvalues of the Hamiltonian. 

There are several possibilities to derive analytically the 
BCS equations from the expectation value (0|iJ pa i r |0) at 



T = 0. One of them is presented, e.g., in Ref. | 
BCS equations are obtained by demanding that 



Ea and c 



0, 



The 



(19) 



where bj and Cj are defined in Eqs. 0. The first expres- 
sion in Eqs. 119fl means that the pairing Hamiltonian is 
diagonalized in the quasiparticles space, the second one 
indicates that the so-called "dangerous diagrams" are ex- 
cluded from the theory. The solution of the BCS equa- 
tions is unique. In another words, it is absolutely neces- 
sary that Eqs. (|19|) are fulfilled exactly to have the BCS 
equations in the form given in Eq. J3J). Let us verify this 
for the MBCS. 

We introduce bj and Cj quantities to replace the 
{uj,Vj} coefficients in Eqs. (0 by {uj,Vj} coefficients 
and calculate the later from the MBCS equations. The 
differences \bj — Ej\ and Cj quantities for several neu- 
tron sub-shells in 120 Sn are presented in Fig. 02a) and 
Fig. HJb) , respectively, showing that 



bj 7^ Ej and Cj ^ 



(20) 



In addition, Eqs. i|19|) are also not fulfilled within the 
MBCS. 

We remind a reader that the MBCS equations were not 
obtained analytically but were written down in analogy 



7 



> 



IW 

i 

IX)" 



0.8 



0.4 



> 

ID 



120 




0.5 1 

T (MeV) 



1.5 



FIG. 8: The MBCS bj and Cj quantities as functions of tem- 
perature for some neutron levels in 120 Sn. See text for details. 



with the T = BCS equations and, as far as the MBCS 
founders had noticed, using a formal similarity in some 
BCS and MBCS expressions. But, as a matter of fact, the 
basic MBCS equations contained in Eqs. (|14|l or Eqs. I|15|) 
cannot be reached from the expectation value (0|if pa ir|0) 
at finite T because Eqs. P0)l. 

One finds in the literature p| that the MBCS equations 
may be obtained from a modified HFB model called the 
MHFB. Here, it was noticed that after application of the 
secondary Bogoliubov transformation of Eq. © , a gener- 
alized particle-density matrix at finite temperature "for- 
mally looks the same as the usual HFB approximation at 
T = 0" 5j. Then, "following the rest of the derivation as 
for the zero-temperature case" p|, the MHFB equations 
were written down. Below we verify thermodynamical 
consistency of the MHFB (MBCS) since Ref. f|] lacks 
such an analysis. 

The total energy of the system in the statistical ap- 
proach has the form 



< H >= Tt(HD), 



(21) 



where D is a density operator and < . . . > means av- 
eraging over the grand canonical ensemble. After grand 
potential minimization, one obtains an expression for the 
system energy as 

£bcs = 2 [C 1 - 2u j) v i + «i] ~ A V G (22) 



in the FT-HFB (or FT-BCS). In the MHFB (or MBCS) 
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FIG. 9: Excitation energy of the PFM calculated as E* (thick 
lines) and < H >* (thin lines) in the MHFB (solid lines) and 
the FT-HFB (dashed lines) with T c = 0.42 MeV for N = 10 
(upper part (a)) and N = 2 (lower part (b)). 



it has the form 



2\",(1 ».,'■ UjVj 



A 2 /G . (23) 



Equation l|23|l appears as Eq. (83) in Ref. [5j or as 
Eq. (45) in Ref. [J. 

In Fig. |3 a) we present the low-T part of Fig. Hfb) (up 
to T c ) with predictions of the MBCS (solid lines) and 
FT-BCS (dashed lines) for the system excitation energy. 
It should be reminded that the results in Fig. 0Jb) were 
obtained for the E* quantity (see Eq. I|16|l) with £ cal- 
culated in the MBCS and FT-BCS from Eq. Q and 
Eq. (|22|l . respectively. These results are plotted by thick 
lines in Fig.[!Ja). The excitation energy calculated in the 
MHFB and FT-HFB as 

< H >* = < H > (T)— < H > (0) 

(in the Gjy = G limit for exact comparison with the 
MBCS and FT-BCS results) is shown by thin lines in the 
same figure. Since -Dmhfb (Eq.(66) in Ref. 0) does not 
equal Dhfb (Eq. (12) in Ref. (5J), the < H >mhfb an d 
< H >h FB quantities are slightly different in Fig. |5| 

In any consistent model the values < H >* and E* 
should not differ significantly since they represent the 
same physical observable. One notices that the accuracy 
of the FT-BCS becomes worse on approaching T c as is 
well known, and disagreement reaches a few percent. The 
picture within the MHFB (MBCS) is completely differ- 
ent: the E^[ BC g quantity is several times smaller than 
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the < H >mhfb quantity almost at any T in this exam- 
ple [13 • We have repeated the calculations in Fig. [?5Ja) 
also for N — 8 and 12 keeping T c fixed. The differences 
between the N — 8, 10, and 12 results are hardly noticed 
by eye for each line in Fig. HJa), i.e. convergence of the 
results in this example is reached with N = 8 — 10. The 
correspondence between < H >* and E* remains still 
acceptable for the FT-HFB (FT-BCS) even for N = 2 
(see Fig. |^b)). On the other hand, disagreement in the 
MHFB (MBCS) results reaches several orders of magni- 
tude in the N — 2 case (notice logarithmic y-scale in 
Fig. Kb)). 

To conclude, thermodynamical inconsistency of the 
MHFB (MBCS) is obvious from this example. This in- 
consistency is detected from T << T c . 

The last result clearly demonstrates that the MBCS 
approach is not justified in the framework of the usual 
statistical approach. This conclusion can be reached by 
other reasoning as well. The point is that the basic 
MBCS equations given by Eqs. (|15fl obtained via the sec- 
ondary Bogoliubov transformation cannot be considered 
as a result of the thermal averaging over the grand canon- 
ical ensemble. Indeed, the bar-quasiparticles in Eq. @ 
and the new correlated ground state 

10} = II + ^<4m^rh) 1°) 

jm 

are temperature dependent although the MBCS founders 
do not use this terminology. A possibility to deduce the 
basic MBCS equations via a variational procedure in ap- 
plication to the avera ge v alue (0|JJ|0) contradicts what 
was proven long ago [l2|: a ground state with such a 
property cannot be constructed in principle in the space 
spanned by eigenvectors of a quantum Hamiltonian (see 
also footnote on page 3 in Ref. 



VI. CONCLUSIONS 

In this article we have studied the modified BCS model 
suggested and explored in a series of papers [S El IS IS • 
We have shown that this model yields many unphysi- 
cal predictions: a negative value for the pairing gap, the 
pairing correlations induced by heating in the closed-shell 
systems, the "superfluid - super-superfluid" phase tran- 
sition. It also predicts that the normal phase does not 
exist at any finite temperature. The MBCS has been 
tested on the picket fence model for which an exact so- 
lution of the pairing problem is available. In addition 
to rather poor description of the exact solutions, it has 
been found that the model severely violates the internal 
particle-hole symmetry of the problem. 

Analysis of the MBCS equations, their derivation, and 
description of different physical observables by this model 
have led us to the following general conclusion: The T- 
range of the MBCS applicability has been determined 
to be far below the conventional critical temperature T c . 
Within this narrow temperature interval the MBCS per- 
formance is worse in comparison with that of the con- 
ventional FT-BCS. Moreover, the MBCS is found to be 
thermodynamically inconsistent. 
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